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Abstract. We prove that, for a smooth two-body potentials, the quantum mean- 
field approximation to the nonlinear Schrodinger equation of the Hartree type is 
stable at the classical limit h — > 0, yielding the classical Vlasov equation. 



Introduction 

Consider a system of iV identical classical particles of unit mass, evolving ac- 
cording to the dynamics generated by the following mean-field Hamiltonian 

N 1 1 

If the two-body interaction <p is sufficiently smooth, the behavior of the system 
for large N is well understood (Seee for instance Ref.s [1],[2],[3]). Namely, if the 
initial positions and velocities of the particles are independently and identically 
distributed according to the probability density / = f(x,v), in the limit iV — > oo 
each particle is distributed at time t according to the distribution / = f(x,v,t), 
independently from the others. Here f(x, v, t) is the solution of the Vlasov equation 



(d t +v-V x +E-V v )f(x,v,t) = 0, 
l 



(2) 



where 

E(x,t) = -V(p*p(x,t) (3) 

and 

p(x,t) = J dvf(x,v,t). (4) 

The quantum analogue of this result has been proved in [4]. The framework is a 
system of N identical bosons interacting by a mean-field potential energy 

U(x 1 ,...,x N ) = —^2<f(xi-Xj). (5) 

Then, if the initial wave function factorizes in the limit N — > oo, each particle 
evolves according to the following nonlinear Schrodinger equation of the Hartree 
type: 

h 2 

(iM t + —A-cp*p)^{t) = (6) 

where 

p(x,t) = \^(x,t)\ 2 . (7) 

and * denotes convolution. Further results concerning the Coulomb interaction (see 
[5], [6]) have been also proved. 

In all these results, however, the convergence is strongly dependent on h. This 
dependence prevents (see the case of the Kac potential below) all situations in which 
N — > oo entails h — > and the system is therefore asymptotically classical. 

In the present paper we address precisely this problem and show that the Vlasov 
equation is indeed recovered in the limit N — > oo even when h — > according to an 
arbitrary law. 

The problem has been approached and solved in Ref. [7] , where the the result is 
obtained via compactness techniques, under the hypothesis that the Fourier trans- 
form of ip is compactly supported. Here we deal with C 2 potentials and compute 
the explicit rate of convergence by means of a constructive method. 
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Our technique is based on the WKB method applied to the iV-particle system. 
In the same spirit of Ref. [8] where the classical limit for a class of nonlinear 
Schrodinger equation has been investigated, we choose the initial wave function in 
such a way that its phase fulfills the classical Hamilton-Jacobi equation. Conse- 
quently the equation for the amplitude, which in our context can be complex valued, 
satisfies an equation of hydrodynamical type. For it we deduce H s estimates via 
the energy method. The regularity estimates we obtain are the key to prove the 
convergence of the Wigner transform of the solution in a rather straightforward 
way. 

1. Problem and results 

Consider an N- particle quantum system described by the following mean-field 
Hamiltonian: 

h 2 N 1 

j=l i<j 

3 d 2 

where Aj = 2 and cp is a two-body smooth potential. 

a =i dx 3\<* 

We are interested in the asymptotic behavior of the system at the limit when 
N — > oo and h — > simultaneously. 

Example: the Kac potential. Consider a system of N identical particles of mass 
m = 1 interacting through the Kac potential 

*<*> = Mx) (L2) 

where A is a large parameter of the same order of iV and (p is a given smooth 
potential. The hamiltonian is: 

h 2 N - 



After the rescaling x = Xq the hamiltonian becomes: 

2 N 



1 / h \ . . _ 

Hn = -2\\) E A ;+ A-1 I>te-^)- ( L4 ) 



^ ft 

Setting \ = N,h= — = — we finally get: 



h 2 N 1 



2 ^ J N 

j=i i<j 



The initial condition ^tv = i/;® iV is assumed to be factorized (and hence sym- 
metric in the exchange of particles). 

For the one-particle wave function we choose initially a WKB state: 

%j){x) = aix)^^ 1 (1.6) 

where a and cr are smooth and independent of h. a is real (more generally, one 
could also consider the case where a ~ J2j>o h j aj(x) is a semiclassical symbol). 
Setting Vl-Ov = ^^<f(xi — £j) > we denote ^jv(-,t) the solution of the Schro- 

dinger equation 

h 2 

{ihd t + —A N - W N )* N (t) = (1.7) 
and introduce its Wigner transform: 

f N (X N ,V N ,t) := 

(1 \ f h h 

— j J dY N e- iYN - VN V N (X N + -Y N ,t)V N (X N - -Y N ,t). (1.8) 

Hereafter we use the shorthand notation = (x\, . . . , xn), Yn = (yi, . . . , un), 
V n = (vi,...,v N )- 



We also introduce the j-particle Wigner functions denned by: 



f]<(X jt V jt t) = (l^Tr) 3 ^) J dX N _^ I dV N -jf N {X N ,V N ,t) 



(i-) 3 ^ dX N _j J dY je - iY r v iy N (X N + ^Y j ,t)<I> N (X N - ^Yj,t) (1.9) 

where X N -j = (x j+1 , . . . , x N ),V N -j = (v j+1 , . . . , v N ), Yj = (y 1 , . . . , yj) and X N + 
hYj = (xi + hyi, . . . ,xj + hyj, Xj+±, . . . , xn). Note that, at time t = 0, fjj* = f® 3 
where fo is the Wigner transform of ip. Hence 

fo ^ p(x)5(v - u(x)) in P' as /i ^ (1.10) 

where p = \a\ 2 and u = Vcr (p = \ao\ 2 in the case where a ~ ^ 7 >o h 3 a j(x)). 



Remark. We remark that /■ is the Wigner transform of the reduced density ma- 
trices: 

p(Xj,Zj) = J dX N _j^{X j ,X N _j)^{Zj,X N _ j ). 
The asymptotics of f^(t) is described by the following theorem: 



Theorem 1.1. Assume if G C 2 (1R 3 ), a G C 2 (R 3 ), d a ip, d a a uniformly bounded 
for \a\ < 2, and a G C 2 fl H 2 (R 3 ). Let h = h(N) -> as N -> oo. T/ien t/iere 
ecczsfc T > 0, sufficiently small, such that for all j G N: 

f 3 N (t)^fj(t) mP'(R^xR^), tG[0,T] (1.11) 

where fj(t) = f{t)® J and fit) is the unique (weak) solution of the classical Vlasov 
equation: 

{d t + v-V x + E-V v )f(x,v,t) = 0. (1.12) 



Here: 



E(x,t) = -Vcp* p(x,t), p(x,t) = J dvf(x,v,t). (1.13) 
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Moreover, fort e [0,T]: 

f(x,v,t) = p(x,t)8{v-u{x,t)) (1.14) 
where the pair (p, u) fulfills the continuity and the momentum balance equations: 

d t p + div(up) = 0; d t u + u ■ Vu = — V0 * p. (1.15a) 
More precisely, for any test function F e ©(M 3 - 7 x IR 3 - 7 ), one has 

(fj N (t) ~ fj(t), F) = 0(h + N- 1 ) (1.15b) 

for h small enough and N large enough. 

We refer to this particular situation as hydrodynamic for obvious reasons. 
Theorem 1.1 is based on some regularity estimates to be established in the next 
section. We will write the solution of the Schrodinger equation (1.7) under the form 

* N (X Nt t) = A N (X N ,ty £!l1 ^ 1 (1.16) 

where S N is the classical action satisfying the Hamilton- Jacobi equation: 

d t S N + l-\VS N \ 2 + W N = (1.17) 

and, consequently, An is the solution of the "transport" equation: 

d t A N + (P N ■ V)A N + -A N divP N = -h-AA N . (1.18) 

Here P N := VS N satisfies: 



d t P" + (P iV • V)P JV = -WW. (1.19) 
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As we shall show later on, the above representation holds for a short time T > 
which may be chosen uniformly in N. 

We then consider initial data which are not of hydrodynamical type. We restrict 
ourselves to classical data, namely those for which the Wigner transform is positive 
and normalized: 

f N (X N ,V N ) = f® N (X N ,V N ) (1.20) 

where / is a classical probability distribution on the one-particle phase space. 

Consider now a one-particle superposition of particular WKB states, given by 
the following density matrix: 

p(x,y) = Jdw e^- {x - y) a(x;w)d(y,w). (1.21) 

The Wigner transform of the density matrix is: 

f p (x,v) = (- 1 -) 3 f dw [ dze lz < v - w) a(x-^z;w)d(x + ^-z,w) = \a(x,v)\ 2 + 0(h). 
2tv J J 2 2 

(1.22) 

Setting \a(x, v)\ = y/J{x, v) we see that / and f p are asymptotically equivalent (in 
V')) in the limit h — > 0. Therefore we assume as initial condition 

JV 

p N (X N ,Y N ) = Hp(xi, yi ) (1.23) 

i=l 

for the density matrix (1.21) with a(x; w) = f(x, w). We will prove 

Theorem 1.2. Assume <p G C 2 (1R 3 ), a G C 2 (M 3 ), d a tp, d a a uniformly bounded 
for \a\ < 2, a compactly supported in w and a(-,w) G C 2 fl H 2 (R 3 ) for all w. Let 
h = h(N) -> as N -> oo. Then for all j G N and t > 0: 
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(1.24) 



where fj(t) = /(t)® J and f(t) is the unique solution of the classical Vlasov equation: 

(dt + vV x + E-V v )f(x,v,t) = (1.25) 



where 

E(x,t) = -V<p* p(x,t), p(x,t) = J dvf(x,v,t), (1.26) 
with initial datum f(x,v) = \a(x,v)\ 2 . 

Note that here the convergence result holds globally in time. 

2. The classical system and its mean-field properties 
Consider the associated Hamiltonian system: 

X N (t) = V N (t); V N (t)i = -1 ^ V V { Xi {t) - Xj {t)) (2.1) 

where 

X N (t) = ( Xl (t)...x N (t)), V N (t) = ( Vl (t)...v N (t)). 

Denote by X^(t, Xn, Vn), Vjv(^) = Vzv(*> Xn, Vn) the solution of the Cauchy prob- 
lem with initial conditions Xn, Vn- For a given initial datum (Xn, Vn), we consider 
the empirical distribution, that is a one-particle time depending measure, defined 
by: 

1 N 

/j N (dx, dv, t) := — S(x - Xi(t))S(v - v l (t))dxdv. (2.2) 
i=i 

The following facts are well known (see e.g. [1],[2],[3]). 

1) fi N (dx,dv,t) is a weak solution of the Vlasov equation (1.12). Namely, for 
any test function h = h(x, v), setting (fj, N (£), h) = J (i N (dx, dv, t)h(x, v), we have: 

f t (» N (t),h) = (» N (t),v ■ V x h) - (p N (t),V x <p * p N (t) ■ V v h). (2.3) 
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This follows by a direct computation. 

2) Weak solutions to the Vlasov equation are continuous with respect to the 
initial datum in the topology of the weak convergence of the measures. 

In particular 1) and 2) imply that if n N (0) — ► / weakly, then fi N (t) — > /(t)weakly, 
where fit) is the unique solution to the Vlasov equation with initial datum /. If / is 
sufficiently regular then f(t) hinerits such a regularity and the solution is classical. 

3) Let f N (X N , Vn, 0) be an initial symmetric iV-particle distribution and let 
f N (X N ,V N ,t) = f N (X N (X N ,V N ,-t),V N (X N ,V N , -t)) be the solution of the Li- 
ouville equation. Define the j-particle marginals by: 

ff^Vj.t) := J dX N _j J dV N -jf N (X N ,V N ,t). (2.4) 

Then, if 

/f - f® J (2.5) 

in the limit — > oo and in the sense of the weak convergence of the measures, 
where / = f(x,v) is a given 1-particle initial distribution, then 

f 3 N (t)^f® j (t) (2.6) 

weakly, where f(t) solves the Vlasov equation with initial condition /. Property 
(2.6) is called propagation of chaos. 

We now specialize the above results to our hydrodynamical case. We suppose 
that initially: 

f(x, v) = f(x, 0) = p(x)S(v - u(x)) (2.7) 
(in the sequel u = Vcr) and denote: 



&(X N ) :=X N {X N ,P N {X N ),t) 
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(2.8) 



where P N (X N ) := {u{xi)}f =1 . Then, for a given test function Fj e V(R 3 x IR 3 ), 

J ffWFjdXjdVj = 

J dX N dV N f® N (X N ,V N )F J (X N (X N ,V N ,t),V N (X N ,V N ,t)) = 
(by the Liouville theorem) 

= J dx NP ® N {x N )F 3 {<s>\x N )\&{x N y) ^ 

f dX.dY,! ' (X j. \ j. I ) /•',( A',. I I (2.9) 

in the limit N — > oo. Here we are using the notation X J N to indicate the vector 
(xi, . . . , Xj) if Xn = (xi, . . . , xn). On the other hand, if f(t) is the solution of the 
Vlasov equation, for a short time t < T it has the form 

/(x, u, t) = p(x, t)8(v - u(x, t)) (2.10) 

with p and w solution of eq.s (1.15) as follows by a direct computation. Indeed if 
the pair (p, a) solves Eq. (1.15), f(x, v, t) given by (2.10) is a solution of the Vlasov 
equation and its uniqueness entails the assertion. 
Moreover observe that the estimate: 

< Ci < \V Xi &(X N )\ < C 2 < 1 (2.11) 

which holds for a short time t < T, with constants C\ and C2 independent of N 
due to the mean-field nature of the interaction, allows us to invert the mapping 
Xn — > $*(Ajv), so that S N , P N , A N exist for such a time interval. The estimate 
(2.11) easily follows from the analysis developed in Section 4. 

We conclude this section summarizing under the form of a Proposition some 
regularity estimates on the classical flow, established in Section 4 below, which will 
be used in the convergence proof. 
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Proposition 2.1. Assume tp G C 2 (R 3 ), a G C 2 (M 3 ), d a <p, d a a uniformly bounded 
for \a\ < 2, and a G C 2 fl H 2 (R S ). Let t G [0,T] wito T sufficiently small. Then, 
there exists C\ > independent of N such that for all N and all j G {1, . . . , N}, 
one has, 

\\V Xj A N (t)\\ L i <d. (2.12) 
Moreover for r G [0, t], denoting by xj, 7 = 1, 2, 3 the components of Xi 

and 

where the constant C is independent of N , k, and i. 

We note that the estimates (2.13) and (2.14) express the weak dependence of 
the position and momentum of the k-th particle with respect to the position of the 
i-th particle at time 0, as it is expected in a mean-field theory. 

3. Convergence 

We are now in position to prove Theorem 1.1. We first observe that, for a time 
interval for which estimates (2.11), (2.12), (2.13) and (2.14) hold, we have classical 
solutions of eq.s (1.17), (1.19) and (1.18). Therefore we can express the j- particle 
Wigner function in terms of A N , S N and P N . For Fj G V(R Sj x R 3 - 7 ), we have: 

/ F l i.X J .\j)J[^.\ J Aj.h<IX l <IV l = 
(J-^j 3 j dX N j dV j j dY j e- iY ^F j {X j ,V j )^ N {X N ^Y j ,t)^ N {X N -^Y j ,t) 

11 



, )t , dX N / dYjFj (Xj , Yj)A (Xn + -Y J: t)A" (X N - -Y^t) 



e ±[S N (X N + ^Y J ,t)-S N (X N -^Y J ,t)]^ 

where Fj is the Fourier transform of Fj in the second variable. Changing variable 
Xn —>■ Xn — ^Yj and using the fact that 

/ dYjsaplFjiXj^WYjl + I dYjSv^\V x .Fj{Xj,Yj)\ < Cj, (3.2) 

J Xj J Xj 

we obtain that (setting Xj = X J N ): 

(3.i) = (^)¥/^/^, W ,y j )^(^ + ^, t) ^(^, t ) 

e i[s N (x N +hYj,t)-s N (x N ,t)} + Q(h). (3.3) 
Now Lagrange's theorem yields: 

r h 

A N (X N + hYj,t) = A N (X N ,t) + / d\V Xj A N (X N + \Yj,t) ■ Y 3 

Jo 

Moreover, since 

HVx^IlL = £ ||V^^||| 2 < Cj; \\A N (t)\\ L2 = \\A n (0)\\l 2 = 1 
i=i 

by (3.2) and Proposition 2.1 we get the estimate 

J dX J dX N J dYj\Fj(Xj,Yj)\\V Xj A N (X N + XYj, t)\\A N (X N ,t)\ < 

h\\A N \\ L2 \\V Xj A N \\ L2 [ dYjSaplFjiX^YjWjlKChy/j. (3.4) 

J x 3 

Hence we can conclude that: 

(3.1) = (J-^j 2 jdXsj dY J F J {X J1 Y J )\A N {X N ,t)\ 2 
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e ±[s N (x N +hY j ,t)-s N (x N ,t)] +0 (h). (3.5) 

Note that 0(h) (as well O(j^) later on) depends on j which however is fixed. 
Furthermore: 

f h 

S N (X N + hYj,t) - S N (X N , t) = / d\P N (X N + XYj, ty ■ Yj = 

Jo 

hP N (X N , t) j ■ Yj + 0(h 2 ) (3.6) 

again by Proposition 2.1. 

Here P n (XnY denotes the projection on the j-particle subspace of the vector 
P N (X N ). 

Hence 

3j 

(3.1) = (^) 2 / dX N J dY J F J (X J ,Y J )\A N (X N ,t)\^ pN ( x -' t y^+0(h) = 

J dX N F(X J ,P N (X N ,ty)\A N (X N ,t)\ 2 + 0(h). (3.7) 
Setting T N = \A N \ 2 we have by (1.18): 

d t T N + div (P N T N ) = B N (3.8) 

where 

B N = l -h(A N AA N - A N AA N ). (3.9) 
The solution of eq. (3.8) has the representation: 

T N (X N , t) = p® Ar ($- t (X 7V )) J N (X N , t) + f dsB N (^- s \X N ))J N (X N , t-s), 

Jo 

(3.10) 

where 

J N (X N ,t)=det\ d *"l XN)) \. (3.11) 

OA N 
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Therefore 

J dx N F j (x j ,p N (x N ,ty)r N (x N ,t) = J dXNp^ix^F^ix^^^x^) 

+ J Q dsJ dX N B N ($-( t - s \X N ),s)J{X N ,t-s)F j {X j ,P N (X N ,t) j ). (3.12) 
The last term in the r.h.s. of (3.12) can be rewritten as: 

J* da J dX N B N {X N ,s)F J {^ t - s \X N )\P N {^ t - s \X N )\t)) = 
y J*ds J dX N {A N AA N - A N AA N ){X N ,s) 

F J (^ t - s \x N y,p N (^ t - s \x N y,t)) 

= Y J ds J dX ^ ANvAN • VF o(- ■■)- A N VA N ■ VFj{. . . )), (3.13) 

here we have integrated by parts and made use of a crucial cancellation. 
We now observe that 

N 

Y,\^„F J (^ t - s Hx N y,p N (^ t - s Hx N y,t))\ (3.14) 
i=i 

can be bounded by a constant dependent on j but not on N. Indeed: 

d x? F J (& t - s \x N y,p N (& t - s \x N y,t)) = 

k=l 1 
Vv k -tj{<P y J {X N )j, Vj,t))\ Vj= pN {<S>{ t-s) {XN)j) — . 



I 



Then by Proposition 1.1 we have that: 



d x? F J (& t - s \x N ) J ,p N (& t - s \x N y,t)) = o(i) 
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if % = 1 . . .j, while 

d x? F(^- s \x N ) j ,p N (^- s \x N ) j ,t)) = o(^) 

if % > j. Hence, by (3.14): 

r - N 
| / dX N A N VA N ■ VFj(. . . )| < \\A\\ L2 |Vx^i(. • ■)\\\^ Xi Ah 2 < Cj. 

J i=i 



Therefore 



(3.1) = J dX N p^ N (X N )F j (^ t (X N ) j , &(X N )j) + 0{h) + 0(1). (3.15) 

Notice that the first term in the r.h.s. of (3.15) is purely classical so that we can 
apply the convergence result (2.9) to conclude the proof. 

The proof of Theorem 1.2 follows along the same lines. Proceeding as as above 
the Wigner function is in this case: 

f N (X N ,V N ) = (I) 3 / dn N j dY N e- iY "- v " 

A N (X N — — Y N , 0, N , t)A N (X N + — Y N , CL N , t) 

e ±[s N (x n -±y n ,q n ,t)-s N ((x n +^y n ,n N ,t)] (3 16 ) 

where A N and S N are the amplitude and the action parametrized by the initial 
momenta Oat. S n and A N are the solution of eq.s (1.17) and (1.18) with initial 
conditions S N (X N , tt N )=tt N -X N and A N (X N , Q N ) = a® N (X N , Q N ). Therefore 
for a test function Fj we have: 



2tt J 



1-1 



) / r/O.v / <IX x I dYjFjiXj.Yj) 
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A N (x N + ^ Yj , n N , t)A N (x N - \y h n N , t^iis"^^^^"^-^^,^. 

2 2 

(3.17) 

Proceeding as in the proof of Theorem 1.1 we find: 

= J dflN J dX]yf® N (Xn, VN)F(XN(t, Xn, Oat), X]y(t, Xn, Oat)) 

+0(h) + 0(±). (3.18) 

Note that that a is assumed compactly supported in w to avoid complications with 
the integral in the initial momenta. However a sufficiently rapid decay of a(-,w), 
for large w, would yield the same result. Of course, once more, everything holds for 
a small time interval. However now the smallness of the time interval depends only 
on the smoothness of the potential <p (because of the particular form a(x) = w ■ x). 
Hence the convergence at time T allows us to extend the argument up to 2T, using 
fact 2) of Section 2 and that f(x,v,T) is still compactly supported in velocity. 
Therefore the convergence can be extended to arbitrary times and the proof of 
Theorem 1.2 is now completed. 

4. Regularity estimates 
In this section we prove Proposition 2.1. We start by considering 

d t P N + (P N ■ V)P N = -WW. (4.1) 

which is independent of 

d t A N + (P N ■ V)A N + -A N divP N = -h-AA N . (4.2) 
to be considered later on. 
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Notice that, for s E [0,t], 



P l N (^ t - s \X N ),t) = ^ t (Y N (s)) l , where Y N (s) = *-(X N ), i = l...N. 

(4.3) 

Using now the short-hand notation Xi(t) = $*(Xjv)i, Pi(t) = ^(X^i and de- 
noting xf(t) and pf(t) the ct-th components, a = 1, 2, 3, we have: 

=Xi+ Pi(s)ds 
Jo 

r l l 

Pi {t) = Va( Xi ) - — VVv?(xi(s) -x fc (s))ds. (4.4) 

where Xtv = (x\ . . .xn)- 

Introducing the force F a = —d x a(p, we have: 



a*2 



= 8i,j8p,7 + ds 



Aft («) 

a*? 



(4.5) 



dpf{t) d 2 a 



dx] " dx]dx? 



J u-ii r* 



dxf(s) dx^(s) 



dx] 



dx] 



(4.6) 



Hence: 

dp?(t) d 2 a 



dx] dxldx? 



nt 1 

kj+ / ds-^^2^2d x c.F^(x i (8)-Xk(8))(Si, j -S k , j )S a ^+ 

J ° k^i a 

dpt(r) dp%{t) 



Jo JO iV !, „ 



(4.7) 



We now observe that, if i ^ j, the first two terms in the r.h.s. of (4.7) are O(jj) 
and hence (taking t G [0, T] with T small enough): 



dx] 



+ 



dpf(t) 



dx] 



(4.8) 
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with C independent of N. 

Remark. Higher derivatives could be handled in the same way to obtain: 



dx] 1 . . . dx] s 

Jl Js 



+ 



d s P Ut) 



dx] 1 . . . dx] s 

Jl Js 



< C 



1 

N 



+ n 6 *>* 



(4.9) 



r=l 



assuming a stronger regularity. 

Furthermore, setting Yjv(s) = (yi(s), . . . , Vn(s)), 

dP?(Y N ,t) 



dx] 



EE 

k f3 



d**(Y N )f , dyl{s) 

Y N =Y N {s)- 



dx] 



(4.10) 



We note that the terms in the sum with k ^ i and k ^ j are O(j^), while the 
term with k = i or k = j are Therefore, if % ^ j, the full sum is O(j^)- 



If % = j the sum is 0(1) because of the term k 
Summarizing: 



0ff ($(*-«> (Xjv),*) 



dx] 



< C 



N 



i = j which is indeed O(l). 
+ 6iJ). (4.11) 



Remark. A similar analysis on the higher derivatives yields: 



d*pf ($(*-«) (X N ),t) 



dx 1 , 1 . . . dx!, s 

Jl Js 



< C 



1 

N 



I'tJr 



(4.12) 



r=l 



We now proceed to analyze eq. (4.2) to obtain estimate of the solution in H 1 . 
Applying the operator V Xj to the equation, we obtain: 

1. 



dtV Xj A + (V Xj P iV • V)i4" + (P iV • V)V K .^ V + -V^divP 
1 



N 



: A N divV Xi P N -h l -AV Xi A N . 



(4.14) 



In computing 
d 



^-(V x .A»,V Xj A») = (d t V Xj A»,V Xj A») + {V Xi A»W Xi A») (4.15) 
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we realize that, due to the symmetry of A, the last term does not give any contri- 
bution. Also, the sum of the terms non involving V X .P N vanishes: 

(V,,^, P N ■ VV Xj A N ) + (P N ■ VV Xj A N , V X] A N )+ 

±(V Xj A N ,V Xj A N divP N ) + ±(divP N V Xj A N ,V x .A N ) = 0. (4.16) 
Here we use the reality of P N and the identity: 

J P N -V\V Xj A N \ 2 = -J divP N \V Xj A N \ 2 (4.17) 

We finally observe that, by eq. (4.11), 

{V X .P N ■ V)A N = (V Xj P 3 N ■ V X .)A N + OiN- 1 \\V Xk A N \\ L2 ) (4.18) 

and thus, denoting WA"^ := (£ fc || V Xk A N \\l 2 )^ (so that £ fc/j l|V, fe ^|| L2 < 
V / iV||^4 JV || i), we arrive to the inequality: 

J^||V^.4.^(*)||i 2 < C-CllV^^^WIli^ -K ^V- 1/2 ||^4. 7Nr || 1 ||V^-4.' Nr (*)IU 2 } (4.19) 

with C independent of N and j. In particular, taking the sum over all j, we obtain, 

±\\A N (t)f 1 <2C\\A N (t)\\l (4.20) 

Since at time zero H^-^ll? is O(N) the same conclusion holds on any time interval. 
Going back to Eq. (4.19), we obtain 

j t \\V Xj A N (t)\\i 2 <C\\V Xj A N (t)\\i 2 (4.21) 

with a new constant C independent of N and j. Then the result follow by observing 

that || V Xj ^^(0)111,2 is uniformly bounded. 
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